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Abstract 

It is pointed out that, contrary to some claims in the literature, the domain walls cannot 
be a source of a correlated at large scales primordial magnetic field, even if the fermionic 
modes bound on the wall had ferromagnetic properties. In a particular model with massive 
(2+1) dimensional fermions bound to a domain wall, previously claimed to exhibit a ferro- 
magnetic behavior, it is explicitly shown that the fermionic system in fact has properties of 
a normal diamagnetic with the susceptibility vanishing at high temperature. 



The existence of magnetic field correlated at a galactic scale is believed to require 
a strong primordial field correlated at cosmological distances at some stage in the early 
universe. This phenomenon would find a natural explanation, if the primordial field was 
created by extended objects, having a cosmological size. One class of extended objects 
that might have existed in the early universe is provided by domain walls, assuming that 
a model containing these walls successfully avoids general constraints on undesirable 
cosmological consequences of such extended defects. There has been several claims in the 
literature P, ^, P] that the modes of fermion field bound to a domain wall produce, in 
certain models, a magnetic field B, which could provide a much needed explanation of a 
primordial magnetic field in the early universe. The purpose of this note is to point out that, 
as enticing as this explanation could be, it cannot be physically correct. Namely, it is almost 
trivial to show that, irrespective of the details of the dynamics, a correlated over the entire 
area of a flat domain wall magnetic field should be greatly suppressed by an inverse power of 
a cosmological size, even if the wall exhibited a ferromagnetic behavior. Furthermore, in the 
models that are sufficiently well formulated the walls are in fact diamagnetic, contrary 

to the previous claims. In this letter the model of Ref. |^ with massive fermionic modes 
with broken parity is considered, and it is shown that the fermion system on the wall is a 
normal diamagnetic, while the previous claim of ferromagnetism is based on an incorrect 
formula p| for the B dependent free energy of the fermionic vacuum. Q 

In order to prove the statement about the suppression of the magnetic field that a domain 
wall might spontaneously create, let us introduce a large normalization box with the sides 
of length L, and consider a flat domain wall, spanning the box along the x, y plane, thus z 
being the coordinate axis perpendicular to the wall. From symmetry, a correlated magnetic 
field can have only the z component B^, which does not depend on the coordinates x and y 
parallel to the wall. Then the Maxwell's equation div B = dictates that B^ also cannot 
depend on the coordinate z, i.e. B^ has to be constant: B^ = B. The energy of such constant 
magnetic field is proportional to the volume of the box, while, the energy, associated with 
dynamics on the wall is proportional to the area of the wall, so that the total energy of 
the system is written as 



with f{B) being the surface energy density in the presence of the magnetic field, determined 

by specific dynamics on the wall. Clearly, for any finite function f{B) the value of B 

iThe claim of || was explicitly shown Q to be based on an incorrect calculation of the energy of the gas 
of massless (2+1) dimensional fermions in magnetic field. 



E{B) = L' — + L'f{B) 
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providing the minimum to the energy goes to zero when the size L of the normahzation 
box is taken to infinity. Thus one concludes that the magnetic field of an infinite fiat domain 
wall has to vanish, independently of the specifics of the model. 

The analyses of the papers |]^, ^, ^ |^ find a ferromagnetic behavior for the function 
f{B) at small field: f{B) — /(O) = —aB, with the spontaneous magnetization a being 
determined by the 'microscopic' parameters of the wall. According to Eq.(|^)the generated 
field then is -B = This consequence of Eq.(0) is acknowledged only in the final version 

of Ref.pi, and is ignored in Refs.@, H, ^. In the latter papers the extent of the field away 
from the wall is assumed to be of the order of the thickness of the wall in direct contradiction 
with the Maxwell's equation. Physically, a finite size L can arise either as the Hubble size 
0, or as a domain size in a multi-domain structure. In either case the magnitude of the 
correlated field is greatly suppressed by and is unlikely to be sufficient for explaining the 
required primordial field. As an example, one can consider the estimates of the field in the 
mechanism of Ref . |p , where the primordial field is associated with the axion domain walls 
at the QCD phase transition, i.e. at the temperature Tqcd ~ ^qcd ~ 0.2 GeV. The size of 
the correlations in the picture presented there is the Hubble scale 1{Tqcd) ~ 30fcm. Thus 
even if one adopts the estimate of Ref. |Q for the spontaneous magnetization: a ~ eAg^^, 
the estimated magnitude of the generated field should be 5 ~ ^ ^qcd/1{Tqcd) ~ lO"^^ 
(instead of the claimed 10^^ G). Furthermore, it will be shown below that the magnetization 
of the axion wall is proportional to the total baryon number density accumulated on the 
wall, so that the estimate of the magnetization is rather on the maximalist side. 

The existence of even a suppressed by magnetic field is contingent on ferromagnetism 

of the fermion modes on the wall. Such behavior was claimed in Refs.0, |, Q within a 

model of massive (2+1) dimensional fermions. Moreover, it has been argued ||^ that the 

ferromagnetic term —B survives in the free energy at all temperatures, in contrast with 

the known behavior in any other systems. In what follows a calculation of the free energy 

of a massive (2+1) dimensional fermion field is presented, in close analogy with a similar 

calculation 0] for a massless case, and it is shown that, contrary to previous claims, the 

system in fact exhibits a normal diamagnetism with the diamagnetic susceptibility vanishing 

at high temperature. The previous erroneous findings of the ferromagnetic behavior were 

due to an inaccurate manipulation with a divergent sum in Ref. [§]. 

^Another simple way to arrive at the same conclusion is to integrate over the surface of the wall the 
magnetic field, created by each element of the surface with constant density of magnetic dipole moment: the 
integral is non-zero only due to edge effects, and vanishes, when the size of the surface goes to infinity. 
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The model of massive (2+1) dimensional fermions is relevant to the situation where the 
fermions have a non-zero mode bound on the wall with eigenvalue m, and there is no mode 
corresponding to the eigenvalue —m. (This un-pairing of the modes arises in a situation 
where the P parity is broken.) The two-dimensional motion of charged fermions in this 
mode is described by a (2+1) dimensional Dirac equation 

ii-f>'id^-teA^)-m) ij = , (2) 

where e is the electric charge, = {Aq, Ai, A2) is the vector potential of an (external) 
electromagnetic field, and 7^ is the set of 2 x 2 gamma matrices, which can be chosen, 
e.g. in terms of the Pauli matrices as 7*^ = a^, 7^ = iai, and 7^ = ia2- It is known (see 
e.g. in Ref.|^) that in an external magnetic field B, perpendicular to the wall, there is an 
asymmetry between positive and negative Landau energy levels, depending on the relative 
sign of eB and m. Namely, assuming for definiteness that m is negative and eB is positive, 
the spectrum of positive energy levels is given by = e B n + im? with n = 1, 2, . . ., 
while that of the negative levels is En- = —\J1e Bn + im? with n = 0, 1, 2, . . .. In other 
words, the level with = is absent from the positive energy part of the spectrum, and is 
present in the negative part. 

The free energy per unit area for the fermion system at a temperature T = l/j3 is 
expressed through the spectrum in the standard way: 

F = F_+F+ + , (3) 

where 

F+ = -p-'— E In (1 + e-"'^™^) (4) 

2^ n=l 

and 

F_ = -/3-'— E In (1 + e-^^2eBnW\ (5) 
2^ n=0 

are the thermal parts of the free energy associated with the real gas of fermions and an- 
tifermions at the Landau levels. The term i^vac in E<l-(0) is the energy of the vacuum state, 
i.e. with all negative energy levels filled and those with positive energy being vacant, 

e B °° 

£^vac = -— E V2e5n + m2 . (6) 
„=o 

The sums for the temperature dependent parts F± are finite, and do not cause controversy. 
The most interesting is the sum in Eq.(^, which is divergent and thus should be handled 
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with some care. In order to make the latter sum tractable without ambiguity, it needs to be 
regularized. A gauge invariant regulator factor should depend on a gauge invariant quantity, 
which naturally can be chosen as the energy of the levels. The exact form of the regulator 
is a matter of convenience, and we use here an exponential form of this factor: exp{—eE^), 
thus making the sum for the regularized vacuum energy read as 

Eil>^{B) = V V2e5rz + m2 exp{-e2e B n - em^) , (7) 

where e is the regulator parameter. In the physically meaningful quantity E^ac{B) — i?vac(0), 
one should take e in the final result. 

The sum in Eq.(^ can be evaluated using Poisson's method based on the identity: 

oo oo oo ^oo 

^ /(n) = / f{x) Y, S{x — n)dx= Y, / f (x) exp{2 IT ikx) dx , (8) 

n=0 •'^ n=-oo k=-oo ^ 

where 5 is an arbitrary number, such that — 1 < 5 < 0. Notice that the sum over n in the 
second expression goes from — oo to +oo. The identity is still valid since the terms with 
< are identically zero. The summand in Eq.(^ is non singular at = 0, so that one can 
in fact set 5 = 0, and write 



eB ^ /-oo 

£; W = V / y/2eBx + m? exp(-e 2 e 5 x - e + 27iikx) dx . (9) 

2^ fc=-oo^O 



In this expression the only term in the sum over k that is singular in the limit e — > is the 
one with k = 0. This term is given by 

/ v2eBx + m? exp(—e2e B x — em^) dx = / Vz + m? exp(— ez — em^) dz 

2tx Jo Att Jo 

(10) 

and does not depend on B. In fact this term is identically equal to the vacuum energy, 
regularized in the same way: 

72. 



^vac(O) = -J ^p' + m2 exp(-ep2 _ ,^2^ Ap_ . (n) 

Thus the k = term in the sum in Eq.@ totally cancels in the difference E^^dB) — -EvaclO), 
and the difference itself is given by the sum of all terms with k ^ 0. The latter sum is finite 
in the limit e — * 0. However an infinitesimal parameter e should be retained in calculation 
of the integrals in order to ensure convergence and proper phase definition of the oscillatory 
integrals. The sum over k ^ is analytically tractable in two limiting cases: large field. 
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eB ^ m^, and small field, eB <^ rn?. In the large field limit the mass can be neglected, and 
the result is given by the massless case p|, 0]: 

E^.,,{B) - E^UO) = (2 e 5)3/2 ^ ^^2) 

io vr^ 

corresponding to diamagnetism with a singular diamagnetic susceptibility. 

The most interesting here is the case of small field, since this is the limit where the 
difference E^s^{B) — _E'vac(0) is claimed H, || to have a linear dependence on B. Using 
the representation in Eq.(P) we find however quite different result. Namely, using Taylor 
expansion in {2eB/m^) for the square root in Eq.(^), and grouping together terms with +/c 
and — fc, one finds in the limit e — 00: 



E^^^{B)-E^^M = -'-^T.T.Tui^^ {-^) ■ 

(13) 

One can easily see that in the sum over s only the terms with odd s are non vanishing, 
thus leaving only the even powers of B in the expansion. Denoting in the non-zero terms 
s = 2p + 1, one finally finds the asymptotic expansion in powers of the field: 

E (B) E rm- ^"B" - (-l)^C(2p + 2)r(|) /ei?y^_ e^B^ 

The leading at small B term in the expansion (0) is positive and quadratic in B, thus 
corresponding to a normal diamagnetism. It is satisfying to verify that, in compliance 
with the common physical intuition, the diamagnetism vanishes at high temperature. By 
"high" one should understand a temperature constrained by the condition T ^ m, since at 
lower temperatures the gas of real fermions cannot produce substantial effect because of the 
thermal blocking factor e"^'*"'. Notice, however that the temperature has to be still lower 
than the mass gap for the fermions on the wall, since otherwise the problem would not be 
reduced to dynamics of a (2+1) dimensional fermion system, but rather one would have to 
consider the full (3+1) dimensional problem. In the latter situation however the effects of 
the wall are small, and no significant phenomena are to be expected. 

In order to calculate the thermal effects under the condition T ^ m ^ e B one can apply 
the Poisson summation formula (|^) for the sums and E_ in Eqs.(^) and (|^). Writing 

^Rearranging the sums is justified as long as the series is absolutely convergent, i.e. at finite e. The limit 
e is taken after the rearrangement. 
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separately the term with k = and grouping together the terms with symmetric non-zero 
values of k results in the expression 

„ „ eB 



TT (3 J 8 

eB 



/ In 1 + exp(-/3 V2eBx + m^) dx - (15) 



{1 OO „QQ 

-ln(l + e-^H)+g / In [i + g-^^^eS^w] ^^2.ik. ^ ^-2.ik.\^ 
2 k=i •'^^ 



Here the lower limit of integration, is such that < 5i < 1, corresponding to summation 
over the Landau levels from n = 1 to oo. (The contribution of the n = level in the negative 
energy spectrum is the first term in parenthesis.) The first term on the right hand side of 
Eq. ( p!5l) arises from the /c = harmonic in the Poisson formula. Taking in this term the limit 
5i +0, one sees that this term does not depend on B and describes the thermal part of 
the free energy of the free fermion gas at zero field: 

d^p 



Fo{f3, m) = -2(3-^ j\n 



1 + exp —[3 Jp^ + rri 



i2ny 



(16) 



The expression in the parenthesis in Eq.(p!5|) vanishes at 5 = 0, since the first term cancels 
against the sum due to the identity 

J2 (e'"''=" + e-'"''=") = -- . (17) 
k=i •'^^ 2 

The rest of the terms can be found by Taylor expansion in B of the integrand in the paren- 
thesis. For the first non vanishing term one finds 



k=l 



2p-P\m\ p2 752 

+0{B^) = - „ , , + OiB^) . (18) 

^ ^ 1 + e-/3IH 24 7r|m| ^ ^ ^ ^ 

Combining this result with Eq. (p!^) , one readily sees that the diamagnetic susceptibility 
indeed vanishes at /3|m| <^ 1, i.e. at T ^ |m|. 

The presented calculation of the free energy of the fermionic system leads us to the 
conclusion that in the considered model the fermionic mode bound to the domain wall gives 
rise to a normal diamagnetism. As usually, the diamagnetic behavior for a system of charged 
fermions arises because the diamagnetism of the 'orbital' motion of the electric charges in the 
magnetic field overcomes the paramagnetism, associated with the normal magnetic moment. 
In order to avoid the diamagnetism, it was suggested [0] that neutral fermions with an 
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anomalous magnetic or electric dipole moment, bound to a domain wall, are fully polarized 
and provide a spontaneous magnetization. Specifically, the strongest effect of the model of 
Ref. arises from a mode for neutrons, bound to an axion domain wall right after the QCD 
chiral transition. Here we present some remarks on this model. 

In terms of a (2+1) dynamics of the bound fermions the Hamiltonian for the anomalous 
interaction has the form 

H^, = ~Ki,a^' F^,i, , (19) 

where k is the anomalous magnetic moment, and a^^ = 1(7^7'' — 'y'^'y^). For a magnetic 
field B in the z direction, B = F12, one has cr^^ = 7", so that the Hamiltonian takes the 
form i^an = —uBip^ip. Thus the spontaneous magnetization a = —{dF/dB)\B=o coincides, 
up to the factor k with the density of the fermionic charge: 

a = K , (20) 

where the averaging over the appropriate thermal state is implied. For the model of Ref . |Q 
the net surface density of neutrons vb = (V^^"^), occupying the zero mode on the axion wall, 
is limited by at least two factors. One is the maximal occupation density for the mode set by 
the requirement that the Fermi energy does not exceed the energy gap A between the modes: 
ub ^ A"^. The value of the gap is model dependent, however, one can perhaps take A ~ Aqcd 
for an estimate. The other limit arises from a consideration of the diffusion of netrons from 
the bulk to the wall. Because of low density of the baryon number in the bulk ub it takes 
time for the netrons to be accumulated on the wall. Taking the mean free path for the 
neutrons at Tqcd ~ ^qcd as ~ 1/Aqcd one estimates the density of neutrons accumulated 
by the wall over the time t as ub ^ tib ^JT/Aqcd- Under a maximalist assumption, that 
the available diffusion time t is given by the Hubble time, t ~ 1{Tqcd), one finds that at 
ub/T'^ ~ 10~^° this limit for vb numerically is also close to Ag^^. 

Summarizing the discussion of this paper, we conclude that the gas of modes of charged 
Dirac fermions on a domain wall is always diamagnetic with susceptibility vanishing at high 
temperature. Therefore such systems cannot spontaneously produce any correlated magnetic 
field. The spontaneous magnetization of the modes for neutral fermions with anomalous 
magnetic moment is proportional to the surface density of the fermion number on the wall, 
and the magnetic field generated by a spontaneous magnetization of a wall is suppressed by 
the inverse of the wall's dimension. These considerations lead us to the general conclusion 
that it is at least highly unlikely that domain walls can provide a physically viable source of 
the primordial magnetic field. 
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